A Detection of Gas Associated with the M 31 Stellar Stream

arXiv: 1506.00635vl [astro-ph.GA] 1 Jun 2015

Accepted for publication in the Astrophysical Journal; May 28, 2015

Preprint typeset using style emulateapj v. 5/2/11

A DETECTION OE GAS ASSOCIATED WITH THE M 31 STELLAR STREAM^

Andreas Koch^, Charles W. Danforth^, R. Michael Rich®, Rodrigo Ibata'^, Brian A. Keeney®
Accepted for publication in the Astrophysical Journal; May 28, 2015

ABSTRACT

Detailed studies of stellar populations in the halos of the Milky Way and the Andromeda (M31)
galaxies have shown increasing numbers of tidal streams and dwarf galaxies, attesting to a complicated
and on-going process of hierarchical structure formation. The most prominent feature in the halo of
M31 is the Giant Stellar Stream, a structure in extent along the sky, which is close to, but

not coincident with the galaxy’s minor axis. The stars that make up this stream are kinematically
and chemically distinct from the other stars in the halo. Here, we present HST/COS high-resolution
ultraviolet absorption spectra of three Active Galactic Nuclei sight lines which probe the M3I halo,
including one that samples gas in the main southwestern portion of the Giant Stream. We see two
clear absorption components in many metal species at velocities typical of the MSI halo and a third,
blue-shifted component which arises in the stream. Photoionization modeling of the column density
ratios in the different components shows gas in an ionization state typical of that seen in other galaxy
halo environments and suggests solar to slightly super-solar metallicity, consistent with previous
findings from stellar spectroscopy.

Subject headings: galaxies: halos — quasars: absorption lines — galaxies: individual (MSI) — galax¬
ies: ISM — galaxies: kinematics and dynamics

1. INTRODUCTION

The assembly of galaxies via hierarchical accretion -
wherein smaller structures merge to create larger struc¬
tures - is now a well-established doctrine of galaxy
formation, which only requires the action of gravity.
This universal prediction of co smological simulations
(e.g., iBovlan-Kolchin et al.ll20I^ is also well observed
in the local universe, with simulations predicting a mul-
titude of substructure s in the Milky Way (MW) halo
(|Helmi fc White! [19991 1. Streams of stars and gas have
been observed in the MW and have b een linked to
previ o us accretion of s a tellite galaxies (iSearle fc Zinnl
119781 : llbata et ^ 11994 iBelokurov et al.l 1200711 . The
Magellanic Stream ( Mathewson et al.l 11974) observed
in 21-cm emission is the result of the interaction be¬
tween the Large and Small Magellanic Clouds (LMC
and SMC, respectively) and the MW and may, gi-
gayears hence, result in the merger of the Magellanic
and MW galaxies. On the s tellar side, the a dvent of
the Sloan Digital Sky Survey (|York et al.l 1200011 has led
to an explosion in the discovery of s t ellar streams in
the MW halo (e.g.. lYannv et al.l[20031 IBelokurov et ahl

2007HGrillmair 1 120091: iNewberg. Yannv. fc WillettI 120(1

i

Martin et al.ll2013h . all of which presumably have their

origin in minor mergers in the Galaxy’s past.

t Based on observations made with the NASA/ESA Hubble
Space Telescope, obtained from the data archive at the Space
Telescope Science Institute. STScI is operated by the Associa¬
tion of Universities for Research in Astronomy, Inc. under NASA
contract NAS 5-26555.

^ Zentrum fiir Astronomic der Universitat Heidelberg, Lan-
dessternwarte, Konigstuhl 12, 69117 Heidelberg, Germany
^ CASA, Department of Astrophysical and Planetary Sciences,
University of Colorado, 389-UCB, Boulder, CO 80309

® Physics and Astronomy Building, 430 Portola Plaza, Box
951547, Department of Physics and Astronomy University of
California, Los Angeles, CA 90095

^ Observatoire de Strasbourg, 11, rue de I’Universite, F-67000,
Strasbourg, France

In the past decade our knowledge of the structure of
our nearest large galactic neighbor, M3I, has under¬
gone a similar tran s forma tion. Deep, wide-held imag¬
ing by llbata et ^ (|2001ll revealed a large stream of
stars extending from near the center of M3I several
degrees to the south east close to the minor axis of
the galaxy. Subsequent studies ((Eerguson et al.l 120021 :
iMcConnachie et al.ll2003l l found that this Giant Stellar
Stream (GSS) exists as a coherent structure more than
4° on the sky, inclined at ^60° with respect to the line of
sight, and extending over > 140 kpc. The GSS is aligned
with the major M31 satellites M32 and NGC205 and,
while it is tempting to associate the stream with those
companions, modeling shows that it most likely repre¬
sents the trailing arm left from a recent, radial merger
of a galaxy of 3 — 5 x 10® M 0 (Mori & Rich 2008, Eardal
et al. 2013). Alternatively it may be the fossil remnant
of a cannibalized satellite or a tidal bridge connecting to
it s neighboring spiral M 33.

iBrown et al.l ()2006l l found from Hubble Space Tele¬
scope (HST) color-magnitude diagrams that the Stream
and halo helds are indistinguishable, which suggests that
only one progenitor was responsible for the debris de¬
posited in inner halo regions. However, if dynamical
mixing were efficient, it could erase the signatures from
different infall sources. Further deep imaging and stel¬
lar spectroscopic studies have shown that the halo of the
M 31-M 33 system is far more complicated than originally
thought and hint s at an eventful formation history (e.g.,
llbata et ^1120071 : IMcConnachie et al.ll2009l l. In addition
to the GSS, many other structures have been discovered
within the halo of M31. Some of these structures are
associated with the large a nd growing population o f An¬
dromeda satellite galaxies (iRichar dson et ^ 1201111 and
globular clusters (|Mackev et al.l 1201011 . while others are
less clearly connected. It has been suggested that many
of the stream-like features in the halo could result from

2

A. Koch et al.

the heating of M31’s thin disk that has been ejected dur¬
ing a recent merger (e.g., Richardson et al. 2008; Bernard
et al. 2015).

We know relatively little about the progenitor of the
GSS. Modeling of the encounter leads to mass estimates
ranging from 3 x 10® Mq (Hammer et al. 2010) to
< 5 X 10® Mq (Mori & Rich 2008); the most recent esti¬
mate of 3 X 10® Mq comes from Fardal et al. (2013), plac¬
ing the system at the mass of the LMC. It has been sug¬
gested that the lack of a compact remnant favors a disk¬
like progenitor. Furthermore, Fardal et al. (2008) argue
that asymmetries in the surface brightness distribution
across the GSS and its morphology in terms of a more
compact core with an extended envelope can best be ex¬
plained if this progenitor disk galaxy was rotating. Based
on the metal-rich nature and the age-metallicity relation
of several halo fields with stream-like color magnitude
diagrams Bernard et al. (2015) suggest an early-type
galaxy as a candidate, which experienced rapid chemical
enrichment to Solar abundances already hy z = 1.

However, no neutral gas analogous to the Magellanic
Stream is evident in the immediate vicinity of the GSS
(Braun et al. 2003). A recent re-analysis of the Hi
data by Lewis et al. (2013) showed the presence of gas
filaments, bridging M31 and the near-by spiral M33,
although they are significantly offset from the stellar
stream. Lewis et al. (2013) also identified a number of
individual gas clouds, but noted a remarkable lack of any
spatial correlation between stellar and gaseous substruc¬
tures, arguing for several dynamical processes at play
that affect both components in different ways. The only
gas associated with the GSS thus appears at the point
where it passes the disk of M31.

Absorption-line spectroscopy provides a different and
complementary tool for analyzing the structure and kine¬
matics of circumgalactic matter (e.g., Rao et al. 2013).
Stellar spectroscopy provides a representative view of
stellar kinematics and overall metallicities, while absorp¬
tion spectroscopy of background sources gives kinematic
and chemical information about the gas phase in the sys¬
tem. In this paper, we present far-ultraviolet (FUV)
spectra of three active galactic nuclei (AGN) obtained
with the Gosmic Origins Spectrograph (COS) onboard
the HST. These three sight lines probe the M31 halo at
projected distances of 26-51 kpc in three different quad¬
rants of the system (Fig. 1).

We see strong absorption from gas near the systemic
velocity of M 31 in all three sight lines. The sight line that
probes the GSS additionally shows clear absorption at
more negative velocities, near the systemic velocity of the
stellar members of the Giant Stream. While Lehner et
al. (2015) associate such absorbers with the hot, gaseous
corona of M31, no connection to the stellar component
of the GSS has been made.

Our measurements of the M31 and GSS-velocity gas
are presented in §3. In §4, we compare the observed gas-
phase detections with the stellar kinematics in the same
field and a kinematic model from the literature. Here,
we also analyze the absorbers by application of basic pho¬
toionization models. Our conclusions are presented in §5.

2. OBSERVATIONS AND DATA REDUCTION

We base our analysi s on HST/COS (jGreen et al.ll2012l :
lOsterman et al.l|2011l) observations of three AGN sight

(kpc)

50 0 -50

—I—,—I—I—I—I—I—I—1—I—I—

-0.60 < [Fe/H] < 0.00

f (degrees)

Figure 1. The locations of the three AGN used in this study
overlaid on a density map of metal-rich stars (after Ibata et al.
2014). Notice that Q0045-|-3926 lies near the edge of the Giant
Stellar Stream. Here, we have masked out the inner disk of M31
within a radius of 17 kpc. The inner halo out to a radius of 30 kpc
and at an axis ratio of 3:5 is displayed at lower contrast to render
visible the inner halo substructure.

lines, which pass near M31 (Fig. 1). Of the three sight
lines, the one near the GSS (Q 0045-1-3926) is notewor¬
thy as it lies 0.2° (~2.7 kpc) East of the main stellar
body of the GSS. The predominance of stellar stream
features and complexity in all regions of the halo is well
documented (e.g. McGonnachie et al. 2009). Three
Seyfert galaxies were observed in 2009, 6-7 October as
part of Program 11632 (P.I. R.M. Rich). Each observa¬
tion consisted of a single exposure in the GOS/G130M
(1135 < A < 1450 A; 6.00 ksec) grating and one to three
exposures with GOS/G160M (1400 < A < 1795 A). In
addition, Q 0045-1-3926 was observed in late 2010 (P.I. N.
Arav). Observation details are given in Table 1.

The exposures were obtained from the Mikulski
Archive for Space Telescopes. The calibrated, one¬
dimensional spectra for each target were next coadded
with the custom IDL procedures described in detail by
iDanforth et al.l (|2010f) . The coaddition routines auto¬
matically scales exposures taken during different epochs
to take into account source variability and stitch G130M
and G160M data together into a continuous spectrum
over the range of 1150-1750 A. Since the two objects
HS 0033-1-4300 and HS 0058-1-4213 were observed with a
single GOS grating position, each dataset features a pair
of ^ 20 A gaps in the spectral coverage, corresponding
to the gap between the two far-UV detectors.

We ht continua to each of the data sets using a semi-
automated line identification and spline-fitting technique
as follows. First, the spectra are split into 5-10 A seg¬
ments. Gontinuum pixels within each segment are iden¬
tified as those for which the signal-to-noise (S/N) ratio

(kpc)

Gas in the M 31 Stellar Stream

3

Table 1

Summary of COS Observations

Target Name

Alternate Name

a

(J2000.0)

5

(J2000.0)

2AGN

Pm 31
[kpc] ^

Date of
of observation

Spectral

Grating

fexp

[s]

Q 0045+3926

lO And

00:48:19.0

+39:41:12

0.134

26.1

2009/10/06

2010/09/07

2010/10/19

GlOOM^’

GieOM^’

GIOOM'^

G130M''

2330

2900

14259

11155

HS 0033+4300

00:36:23.0

+43:16:40

0.120

31.8

2009/10/06

GIOOM^’

GieOM^’

2455

9145

HS 0058+4213

J 0101+4229

01:01:31.2

+42:29:35

0.190

50.7

2009/10/07

GIOOM*"

GieOM*"

2320

2970

^ Minimum impact parameter from sight line to center of M31 based on the scale 228 pc arcmin ^ and mean
distance 792 di 193 kpc.

^ Program 11632, P.l. Rich
^ Program 11686, P.L Arav

vector is less than 1.5 a below the median S/N value
for all the pixels in the segment. Thus, absorption lines
(flux lower than the segment average) are excluded, as
are regions of increased noise (error higher than segment
average). The process is iterated to convergence and a
spline is fitted through the mean continuum pixel value in
each segment. We check the continuum fit of each entire
spectrum manually and the continuum region identifica¬
tions are adjusted as needed. The continuum identifica¬
tion and spline-fitting processes work reasonably well for
smoothly varying data, but they were augmented with
piecewise-continuous Legendre polynomial fits in a few
cases. In particular, spline fits perform poorly in re¬
gions of sharp spectral curvature, such as the Galactic
Lya absorption and at the peaks of cuspy emission lines.
Complete details on the process are given in Danforth et
al. (2014).

The region around M 31 has previously been mapped
in the 21 cm emiss ion line of H I by several authors.
iThilker et al.l (|2004ll conducted a Gree n Bank Telescope
GBT ) survey of the region, while iBraun &: ThilkeJ
2004fl present a deep survey in the M31 region made
with the Westerbork Synthesis Ra dio Telescope (WSR T)
operated in total-power mode (see IBraun et al.ll2003l for
details). While the absorption-line data are much more
sensitive than emission, the COS data give no informa¬
tion on the neutral hydrogen column density since the
MW damped Lya absorption profile extends over thou¬
sands of km s“^.

Unfortunately, the velocity calibration of the UV spec¬
tra is uncertain to a much larger degree than radio emis¬
sion data. HST spectra are recorded in a heliocentric
velocity frame with an on-board wavelength calibration
lamp. However, the precise wavelength solution varies by
^0.2 A (^2 resolution elements) over the detector and
the zero point is uncertain by ^^15 kms“^ due to inac¬
curacies in centering the target in the COS aperture in¬
herent to the target acquisition procedure. The coalition
process cross-correlates strong absorption lines in each
detector segment in each exposure with each other be¬
fore coaddition, so the relative wavelength solution over
a few Angstroms is reasonably accurate.

However, no cross-segment wavelength correlation is
attempted. We correct for these effects to first order
by aligning narrow low-ionization metal absorption lines
in the data in velocity space to the H I 21-cm profiles

from the Leiden-Argentine-Bonn survey (|Kalberla et al.l
[200^ . The H I emission toward all three sight lines
shows a peak at 0 < ulsr < +3 km s“^ with significant,
distributed emission toward negative velocities. Fits to
the Q 0045-1-3926 COS data show the reddest absorp¬
tion component in several ions at ^ —20 km s“^. Poor
data quality in the other two spectra hampers this analy¬
sis rather severely, but the strongest Galactic absorption
appears to be at —40 ^ v < —30 km s“^. These values do
not vary considerably with rest wavelength, so we adopt
a uniform velocity shift of -1-20 km s“^ for Q 0045-1-3926
and -k35 km s”! for HS 0033-^4300 and HS 0058-^4213
to bring the data approximately into the Local Standard
of Rest (LSR) reference frame. Overall, we estimate a
systematic velocity uncertainty of ±15 km s“^ inherent
in the absolute velocity of any absorption measurements
henceforth.

We measure absorption features in three ways.
Firstly, unsaturated absorption features with clear com¬
ponent structure can be htted via the apparent optical
dept h (AOD) and prohle-weighted mean quantity meth -
ods (|Savage fc SeinbadilllOOltTSembach fc Savaggll992ll .

Secondly, multi-component Voigt profile fits often pro¬
vide more accurate measurements in cases of mildly sat¬
urated lines. The significance of any individual line fit
is calculated based on the local S/N ratio, the Doppler
b-parame ter, and fitted equiv alent width via the algo¬
rithms of iKeenev et al.l (|2012ll . Finally, we make use of
simultaneous multi-component Voigt profile fits to mul¬
tiple transitions of the same species, e.g., Si II 1260.42,
1193.29, 1190.42, and 1526.70 A. These techniques pro¬
duce equivalent results for absorption profiles of mod¬
erate strength, but each has its advantages in different
circumstances: Voigt fitting, which is generally preferred,
can accurately measure lines which are moderately sat¬
urated (r ~2-3). The down-side of this technique is
that a velocity component structure must be assumed.
AOD measurements are well-suited for shallow, noisy ab¬
sorbers or where the component structure is uncertain,
while this method tends to under-estimate column den¬
sities for very strong absorbers. The latter does not pose
a concern in our only mildly saturated absorbers. For
details on the consistencies between the methods, we re¬
fer the reader to Danforth & Shull (2008; and references
therein). Measurements from all three techniques are
performed on the data and the resulting values are pre-

4

A. Koch et al.

sented in Tables 2 and 3, where we label the absorber
“components “ by the error-weighted mean of all species.

3. RESULTS
3.1. Q 0045+3926

The spectrum of Q 0045-1-3926 shows five distinct ab¬
sorption components in the range —400 < vur ^
0 km s“^ in many metal absorption lines from low-
ionization Al II and Fe II to highly-ionized C IV tran¬
sitions (Fig. [5]). The strong absorption at n « 0 kms“^

-600 -400 -200 0 200

Velocity (LSR, km/s)

Figure 2. Normalized absorption profiles in various transitions
along the Q 0045+3926 sight line. Data are binned to three pixels
(~40% of the COS resolution) in all cases. The (small) normalized
error vector is shown in gray. Vertical dashes correspond approxi¬
mately to the velocity component centroids discussed in the text.

is similar to that seen in all other extragalactic UV ab¬
sorption lines and can be unambiguously associated with
absorption from the MW disk gas. Other absorption,
often blended with the MW gas, appears at —70
km s“^. This component could plausibly be associated
with the M31 halo given the overlap of MW and M31
stellar velocities discussed in Sect. 4.1 below. However,

this absorption is also consistent with infalling or out¬
flowing gas within the Milky Way halo in the form of
the high velocity clouds (HVCs) and we deem this to be
a more likely explanation (e.g., Wakker & van Woerden
1997; Lehner et al. 2012).

At n = —170 and —230 km s“^, two strong absorp¬
tion components can be seen in transitions of Si II (1193,
1190, 1526 A), Si III (1206.5 A), Si IV (1393, 1402 A),
and C IV (1548, 1550 A). The strong Si II 1260 A line
is blended with weaker S II 1259 A absorption from the
MW at that velocity, though the observed profile is con¬
sistent with the same absorption seen in three weaker
transitions of the same ion. The strong C II 1334.5 A
line is similarly complicated by HVC absorption from
the excited-state C II* 1335.7 A transition, but is again
consistent with the absorption component structure seen
at u ~ —230 km s“^ in other metal ions. Notably, the
relative strengths of the two lines in low ions as com¬
pared with higher ions in the same species (e.g., Si II
and Si IV) show opposite trends implying that the ma¬
terial at n ~ —170 km s“^ is less ionized than the gas at
V —230 km s“^.

Prominent absorption is seen at n ~ —370 km s“^ in

Si II 1260 A, C II, and Si III and weakly in other Si II
transitions as well as both Si IV lines and C IV 1548 A.
This is a strong, distinct, system that is ^ —70 kms“^
relative to the M31 systemic velocity, and we associate
it with the GSS, despite the lack of obvious H I emission
in the immediate GSS, given the lower column densities
that our observations sample. We will further discuss its
possible origin and connections to stellar kinematics in
its surroundings in Sect. 4.1.

We note that all the detected three high velocity com¬
ponents could in principle arise in Milky Way HVGs as
well, although we deem this very unlikely: MW halo
gas clouds are typically observed at velocities up to
^ 200 km s“^, and the ex treme wing of the distribution
reaches up to 400 km s“^ (|Wakker fc van WoerdenifIMI:
iPutman et al.l 1^00211 . Some weak metal-line absorption
systems have been detected in the MW halo that are not
associated with 21-cm emission. These detections have
come in both low-ioniza t ion sp ecies such as Ga II and
Na I (iRichter et al.ll2005l l2009li and higher ions, such as
G IV and O VI ( Sembach et al.l[^003ll .

3.2. HS 0033+4300 and HS 0058+4213

The data for HS0033-f4300 and HS0058-f4213 are
considerably noisier than those of Q 0045-1-3926 due to
fainter source fluxes and shorter exposure times. De¬
tailed hts to the absorption profiles in these two sight
lines are difficult, but similar trends are qualitatively
seen. Fig. [3] shows absorption aX v -1-15 and —70
km s“^ in all bands which we ascribe to MW disk and
halo gas. Both sight lines show absorption at n ^ —180
to —170 km s“^ in the moderately and highly-ionized
species as well, presumably tracing M 31-velocity gas de¬
spite their considerable distance from the main body of
the disk. Both the MW and M31 absorption may plau¬
sibly be sub-divided into additional components as seen
in the higher-quality data of Q 0045-1-3926.

HS 0033-1-4300 shows a low-significance absorption fea¬
ture in several ions (Si III, G IV) consistent with v ^

Gas in the M 31 Stellar Stream

5

Table 2

Species Measurements Toward Q 0045+3926

Species

Velocity

[kms“l]

6-value
[kms“^]

log N
[cm-2]

Notes

V ^ —170 km s ^ component

OI

-170 ±5

20 :

12.7 I*"

< Ifj, AOD, Voigt

Sill

-158 ±4

21 ± 7

13.39 ±0.06

A1190, 1193; AOD, Voigt

Si III

-162 ±2

19 ±4

13.17 ±0.11

AOD, Voigt

Si IV

-160 ±7

22 ±5

12.92 ±0.12

doublet; AOD, Voigt

CII

-130 :

30 ;

14 :

Voigt only, strongly blended

CIV

-173 ±2

18 ±3

13.57 ±0.03

AOD (doublet), Voigt (A1548 only)

NV

< 12.97

doublet; AOD upper limit

V ^ —230 km s ^ component

OI

-236 ± 7

20 :

13.1 I*"

AOD, Voigt la)

Sill

-222 ± 16

24 ±7

12.99 ±0.29

A1190, 1193; AOD, Voigt

Si III

-228 ±4

30 ±5

13.16 ±0.04

AOD, Voigt

Si IV

-225 ±4

28 ±4

13.22 ±0.10

doublet; AOD, Voigt

CII

-227 ±3

24 ± 5

13.89 ±0.02

AOD, Voigt

CIV

-230 ± 3

25 ±7

13.98 ±0.07

doublet; AOD, Voigt

NV

-230 ± 5

10 :

12.97 ±0.12

~ 3(7, doublet; AOD, Voigt

V —370 km s ^ component

OI

-374 ± 7

20 :

13.4 I*"

AOD, Voigt 2a)

Si II

-359 ±4

24 ±2

12.88 ±0.10

A1190, 1193,1260; AOD, Voigt

Si III

-355 ± 2

35 ±6

13.09 ±0.04

AOD, Voigt

Si IV

-390 ± 9

40 :

12.62 ±0.11

A1393 only; AOD, Voigt

CII

-351 ± 5

29 ±6

14.12 ±0.05

AOD, Voigt

CIV

-378 ± 1

34 ±5

13.20 ±0.01

AOD (doublet), Voigt (A1548 only)

NV

-393 ± 16

30 :

12.65 ±0.06

< 3a; A1238 only; AOD, Voigt

^ All line profiles have been shifted by +20 km s ^ as discussed in the text.

^ The 3a detection limit for 01 lines is log N < 13.54. Marginal, lower-significance mea¬
surements are presented for completeness.

Table 3

Species Measurements Toward HS 0033+4300 and HS 0058+4213

Species

Velocity

[kms“^]

6-value
[kms“^]

log N

Notes

[cm J

HS0033±4300, v

^ —170 km s

^ component

OI

< 14.15

upper limit

Sill

< 13.2

A1190,1193,1526; simulfit

Si III

-142 ± 3

68 ±4

13.5 ±0.3

AOD, Voigt

Si IV

-177 ±2

53 ±18

13.5 ± 0.1

A1402 only; AOD, Voigt

CII

-155 :

65 :

14.5 :

very uncertain, Voigt, AOD

CIV

-176 ± 4

35 ± 7

14.13 ±0.07

doublet; AOD, Voigt

NV

< 13.56

upper limit

HS0033±4300, v

~ —350 km s“

^ component

OI

< 14.15

upper limit

Sill

< 13.2

upper limit

Si III

-350 ± 5

22 :

12.9 :

< 3(t Voigt, AOD

Si IV

A1402; instrumental feature

CII

-359 ± 11

> 50

14.3 ±0.1

AOD, Voigt

CIV

-323 ±21

75 ±12

13.77 ±0.09

doublet; AOD, Voigt

NV

< 13.56

upper limit

HS0058±4213, v

^ —180 km s“

^ component

OI

< 13.9

upper limit

Sill

-188 ± 11

37± 15

13.8 ±0.2

A1190, 1193,1526; AOD, Voigt

Si III

-179 ± 10

15 ±5

14.5 ± 0.3

Voigt

Si IV

-183 ± 8

38 ±7

13.45 ±0.18

doublet; AOD, Voigt

CII

-183 ±2

33 ±6

14.5 ± 0.1

AOD, Voigt

CIV

-162 ± 8

53 ±14

14.15 ± 0.13

doublet; AOD, Voigt, simulfit

NV

< 13.4

upper limit

^ All line profiles have been shifted by +35 km s ^ as discussed in the text.

6

A. Koch et al.

-600 -400 -200 0 200 -600 -400 -200 0 200

Velocity (LSR, km/s) Velocity (LSR, km/s)

Figure 3. Same as Fig. [2] for HS 0033+4300 (left) and HS 0058+4213 (right). Data are binned by three pixels. The stronger Si IV line
(A = 1393) is blended in the left panel with strong intrinsic absorption, so the weaker A = 1402 line is shown. Vertical lines are the five
components seen in Q 0045+3926.

Gas in the M 31 Stellar Stream

7

—350 km gas, though there is little evidence for
this absorption in lower-ionization species. This con¬
trasts Q 0045-1-3926 where the strongest absorption at
this velocity is seen in the low-ionization species Si II,
Si III, and C II. No convincing absorption is seen in the
HS 0058-1-4213 sight line at n < —300 km s“^.

4. DISCUSSION

Our examination of the three absorption spectra finds
that all have in common Galactic components, and com¬
ponents at around —170 to —230 kms“^ that are plau¬
sibly associated with the M31 disk.

4.1. Gas Kinematics and Comparison with Stellar
Samples

Q 0045-1-3926 stands out in exhibiting a unique com¬
ponent at —370 kms“^ not seen in the other two spectra,
although there is a low-significance absorption at —350
kms“^ toward HS 0033-1-4300. The simulations of Mori

6 Rich (2008) suggest that a stream component would
be found in this NW region, but the putative components
appear too weak for a firm association. As the sightline
towards Q 0045-1-3926 passes ^2.7 kpc from the GSS
we now turn to the question of whether observed stel¬
lar radial velocities or modeled stellar stream velocities
are consistent with the —370 kms“^ component in the
gas. Fig. 4 (bottom panel) shows stellar radial veloci¬
ties derived from Keck/DEIMOS measurements within

7 kpc of Q0045-1-3926 (Ibata et al. 2015 in prep.); the
stream shows a clear signature with a peak near —500
kms“^, well separated from the gas component. Such a

Velocity v [km s

Figure 4. Observed stellar radial velocity distributions in the
immediate vicinity (7 kpc) of two of our AGN sightlines, where
available: HS 0058-1-4213 (top panel) and Q 0045-1-3926 (bottom
panel), taken from Ibata et al. (2004) and Ibata et al. (2015, in
prep.). The bottom field contains a strong signature of the Giant
Stream at —500 kms“^. Solid and dashed lines indicate the detec¬
tion of gas absorption at those velocities, as listed in Tables 2,3.

difference is in line with the lack of any obvious spatial

correlation between stellar and gas velocities reported
by Lewis et al. (2013). The presence of a strong stellar
velocity gradient within the GSS across its 4.5° extent
has already been noted by Ibata et al. (2004), ranging
from —300 kms“^ near its tip to —500 kms“^ towards
the disk of M31 (see also Koch et al. 2008; Gilbert et al.
2009). While the gas filaments also display a strong kine¬
matic gradient, there are two marked differences (Lewis
et al. 2013): firstly, gas is only detected with a signif¬
icant parallel offset to the stellar stream and secondly,
gas near the disk shows the systemic velocity of M31, at
—300 kms“^, transiting to more positive velocities as it
approaches M 33. Koch et al. (2008) and Gilbert et al.
(2009) also identify a kinematic signature of the GSS at
—440 kms“^ and thus much closer to the gas velocity
of our present work. The respective fields (“a3”) are,
however, located at ^30 kpc from Q 0045-1-3926, adding
further complexity to the interpretation of velocity gra¬
dients and offsets between stream and stellar dynamics
and spatial coincidences.

The top panel of Fig. 5, in turn, shows the Mori & Rich
(2008) model distribution of velocities for the M 31 halo
and stream particles within 5 kpc around Q 0045-1-3926.
The models had been tailored to reproduce the morphol¬
ogy and kinematics of the GSS and distinguished parti¬
cles from the disk/bulge, halo, and the accreted satellite.
First, we note that no disk or bulge particles, ejected
during the merger, are found in this region. Halo stars
in the simulation show, per construction, the broad dis¬
tribution of stellar velocities that also include the three
absorbers we observe in the gas phase and in the stel¬
lar samples. Turning to the satellite particles that make
up the stream in the simulations, there is a clear coin¬
cidence of the broad distribution with the observed gas
velocities at —370 km s“^. Moreover, the stream ve¬
locity distribution shows kurtosit0 toward less negative
velocities with a peak at —300 kms“^.

We now turn to the case of HS 0058-1-4213. Fig. 4 (top)
shows that the respective stellar sample was dominated
by MW foreground dwarfs above ^ —100 kms“^, which
level off into the M31 extended disk and halo profile be¬
low ^ —150 kms“^; a tail of stars is found at —400
kms“^ < V < —300 kms“^, which represents the canon¬
ical halo population. The observed gas absorption in this
sight line at —180 kms“^ is also consistent with veloci¬
ties of individual gas clouds in the M31 halo (e.g., Lewis
et al. 2013) so that both stellar and gas velocities in
this region probe M31’s halo with no obvious connection
to any substructures. We note that no Stream particles
are found in this region of the Mori & Rich (2008) sim¬
ulations and the only dominant component is the broad,
halo model peak at —300 kms“^.

Finally, while no stellar velocities are available in the
vicinity of HS 0033-1-4300, the simulations predict that a
rather diffuse stream component should be seen at highly
negative velocities and another peak that blends into the
MW foreground, above ^0 kms“^ (Fig- 5, bottom). The
gas absorption seen at —170 and, more weakly, —350
kms“^ thus appears to be consistent with typical halo
gas.

® Prom the 140 stream particles within 5 kpc of Q 0045-1-3926,
we find a mean velocity of —340 kms“^ with a Icr-dispersion of 40
kms“^, skewness of 1.2, and kurtosis of 3.1.

A. Koch et al.

0.25
0.2
0.15
o 0.1

.Q

I 0.05

Z

0) U

N

1 0.06

o

z

0.04

0.02

-800 -600 -400 -200 0 200

Velocity v [km s“^]

Figure 5. Same as Fig.[^ but fo r the r adial velocity distributions
from the models of IMorifc Richl ()2QQ8l h The gray-shaded and
clear histograms are drawn from the halo, and stream particles,
respectively, while vertical lines indicate the detection of gas ab¬
sorption at those velocities in the respective sight lines, as listed
in Table 2. No stream particles are present in the vicinity of HS
0058-h4213.

All three sight lines thus show absorption at near
the systemic velocity of the M31 halo (v ^ —170
km s“^), considering the stellar halo’s broad velocity dis¬
persion (at ^85 kms“^; Reitzel & Guhathakurta 2002).
Q 0045-1-3926 shows two clear components separated by
^ 60 km s“^ and it is clear from the ratios of high
ions to low ions that these components have different
ionization properties in different sightlines, in that the
V ~ —170 kms“^ system shows weaker CIV and stronger
low-ionization lines than the one v ^ —230 kms“^. Sim¬
ilarly, towards HS 0033-1-4300, C IV and Si IV absorption
is seen near v ^ —176 km s“^ while C II and Si III show
absorption centroids at —150 km s“^.

In conclusion, we note that stars are found with veloc¬
ities consistent with all the three absorbing systems in
the vicinity of the AGN sight lines, which is a manifes¬
tation of the large velocity dispersion of the M31 halo
system. Based on the observed stell ar velocities and the
comparison with the simulations of IMori fc Richl (|2008ll
we argue that the lowest-velocity component present in
the AGN spectra are associated with the Stream itself.

4.2. Modeling Physical Conditions

With our measurements of many metal ion column
densities, it is possible to constrain the physical condi¬
tions in the absorbing gas with simple photoionization
models. To this end, we compare measured column den¬
sities with a grid of GLOUDY models dev eloped for cir-
cumg alactic absorbers (described fully in iKeenev et al.l
1201311 . Briefly, a slab of gas is illuminated by an ion¬
izing c ontinuum given b y the extragalactic radiation
held of iHaardt fc Madaul (120121) . As in Keeney et al.
(2013), we use a purely-extragalactic (AGN) background
since all three of our sight lines probe regions outside
M31’s proximity zone where ionizing hux from early-

-Stream

-Q0045+3926

-Stream

-HS0033+4300

type stars is expected to have an effect (Uprox 20 kpc;
iGiroux fc Shulll [1997I1 . We calculate a grid of models
where free parameters ionization parameter U = n^/nn
and metallicity Z are varied to provide a grid of model
conditions. Best-Ht model parameters for the three ab¬
sorption components toward Q 0045-1-3926 are listed in
Table 4.

Measurements in three adjacent ionization states of a
metal provide a strong constraint on the ionization pa¬
rameter. In these data the strongest constraints are pro¬
vided by Si II/III/IV, though GII/IV ratios are generally
consistent.

H I absorption at the velocity of M31 is unrecover-
ably blended with Galactic absorption which makes the
metallicity of the absorbing gas difficult to constrain di¬
rectly. We set an upper limit of Nui < 5 x 10^^ cm“^
since the Q 0045-1-3926 sight line lies out side the low¬
est H I 21 cm column density contour of iThilker et al.l
(|2004f) . However, we can infer Vhi for each compo¬
nent by hnding the range of total column density Nh
over which all of the silicon ion column densities in a
given component can be explained simultaneously in a
single, purely-photoionized phase of a given metallicity
Z. This gives values for Nh and Z. These quantities,
along with the hydrogen neutral fraction {fni, uniquely
given by the ionization parameter above) determines the
neutral hydrogen c olumn density Ah t all within fairly
large uncertainties (IStocke et al.lfeOOTl) . The equilibrium
temperature, T, is a value returned by the best-fit U, Z
model. A line-of-sight cloud size is then inferred in turn
as Del = riH Jhi /Nbi ■ Assuming a spherical cloud of
constant density, we calculate a cloud mass Md-

We note that many assumptions have gone into these
models and quantities farther down the chain of logic
should be taken correspondingly with appropriate cau¬
tion. Modeled ionization parameters are probably rea¬
sonably accurate while inferred cloud mass may be inac¬
curate by orders of magnitude or more.

The GLOUDY model best fit parameters are log U k,
—3.0 for both the GSS gas and the component at —170
km s“^, typical of conditions seen in the MW halo (e.g.,
Shull et al. 2011). The gas at —230 km s“^ appears to be
more highly ionized (log U ~ —2.6) and thus typical of
gas in “highly-ionized” HVGs, which is in agreement with
the qualitative estimates of the respective line strengths
seen in Fig. 2.

The derived metallicities for the three absorption com¬
ponents along the Q 0045-1-3926 sight line are all consis¬
tent with solar or slightly super-solar metallicity, though
with large uncertainties. In particular, for the two M 31-
velocity components which are kinematically unrelated
to the GSS and presumably trace gas in the halo of M 31
(see also Lehner et al. 2015), these values are higher
than the average stella r metalliciti e s in t he, generally
more metal-poor, halo. iKoch et ^ (1200811 found stellar
metallicities of [Fe/H]ss —1.4 dex in fields at 20-40 kpc
from the galactic center. Stars in the Giant Stream, how¬
ever, are known to be more metal-rich than the surround¬
ing halo stars (Brown et al. 2006; Koch et al. 2008).
For instance, Ibata et al. (2001) reported on an average
[Fe/H] of approximately —0.7 dex, while Bernard et al.
(2015) list a median metallicity of —0.35 dex based on
their HST photometry. Moreover, the metallicity distri-

Gas in the M 31 Stellar Stream

9

Table 4

Derived Physical Quantities from CLOUDY Models

Absorber

log AhG
[cm-2]

log U

log Z/Zq^

log riH
[cm“^ ]

logT

[K]

log /hi

Del

[kpc], line of sight

log

[Mq]

Q 0045-1-3926, -370 km s”!

15.96 ±0.33

QC + O.IS

-0.13t°;34

-3.02

4.04

-2.06

0.35

2.86

Q 0045-1-3926, -230 km s”!

15.36 ±0.30

o f?o+0.26

±0.48/°;^"

-3.45

3.77

-2.28

0.40

2.58

Q 0045-1-3926, -170 km s”!

16.02 ± 0.48

— 2

24

±0.l8l°;«3°

-3.08

3.83

-1.96

0.38

2.88

HS 0033-1-4300, -350 km s"!

17.3 ±0.4

-3.o1°;3

_1_0 1+0-3

-3.1

4.2

-2.3

17.0

7.8

HS 0033-1-4300, -170 km s”!

15.9 ±0.9

_1_0 2^^'^

-3.6

4.2

-2.7

4.9

5.7

HS 0058-1-4213, -180 km s”!

17.4 ±0.3

-3.1 ±0.1

—0 7+0-3
‘ -0.2

-2.9

4.1

-2.0

8.6

7.1

A^hi is determined by assuming that Sill/III/IV all arise in a single photoionized phase.

For this, CLOUDY adopts the Solar metallicity of log Zq =0.012 from Asplund et al. (2005).

Mci is the cloud mass (in solar masses) assuming a spherical cloud with density riff and diameter D^i-

butions of Gilbert et al. (2009) peak at —0.2 dex and in¬
dicate a significant population of stars with roughly Solar
metallicity in fields dominated by GSS material. These
findings support our measurement of Solar metallicity in
the stream gas, especially if one considers all the dif¬
ferent, employed measurement techniques in the studies,
from stellar calcium triplet spectroscopy, to color magni¬
tude diagram fitting, and our modeling of the gas phase
metallicity.

Given the poor quality of the data for the
HS 0058-1-4213 and HS 0033-1-4300 sight lines we hesitate
to place much weight on their model results. However,
they are presented in Table 4 for completeness.

5. SUMMARY

We have investigated three sight lines in the M31 halo
using HST/GOS and AGN as background sources. All
of the sight lines exhibit components that we associate
with the Milky Way, and as well, a component at —170
to —230 kms“^ that we assign to the disk of M31. In
the sight line of Q 0045-1-3926, which is projected ^2.7
kpc from M31’s Giant Stellar Stream, we find S II and
C II absorption at —370 kms“^, the most negative de¬
tected velocities in our sample. Here, we argue that
the —370 kms“^ component in Q 0045-1-3946 is certainly
associated with structures in M31, most likely with the
GSS.

The bulk of stellar velocities in the GSS covers a broad
range, reflecting a strong radial velocity gradient across
its body. Stars in the immediate vicinity of the AGN
sightline move at near —500 kms“^, while other stream
components have been found at velocities closer to our
observed gas velocity (Ibata et al. 2004, 2015; Koch et al.
2008; Gilbert et al. 2009), albeit at larger projected dis¬
tances. Aided by simulations of the merger that formed
the GSS, Fardal et al. (2008) suggested that many of
its morphological features could best be reproduced with
a rotating disk-progenitor. In that case, a detachment
of the stellar from the gas velocities would occur nat¬
urally. It is obvious that different mechanisms were at
play, contrasting the tidal disruption of the stellar com¬
ponent of the stream progenitor with the gas removal
by ram-pressure stripping in M31’s hot circumgalactic
corona (Mayer et al. 2006, Lewis et al. 2013; Lehner
et al. 2015). Likewise, the remarkable lack of coherence
between stellar structures and presently detected large-
scale gas in the halo emphasizes the complexity of the
GSS formation and dynamics (Lewis et al. 2013).

Further evidence for a stream-gas connection comes
from the putatively high metallicity we derived for the
gas phase. The halo of M 31 has long been known to have
stars with metallicity similar to that of 47 Tuc (Mould
& Kristian 1983; Rich et al. 1996). Even higher metal-
licities were found in larger surveys, allowing for a bet¬
ter characterization of the main halo polluters. In fact,
the HST-based age-metallicity relation of Bernard et al.
(2015) showed that halo fields rich in substructures ex¬
perienced very rapid chemical enrichment so that, by a
redshift of z = I, Solar metallicities were in place. This
kind of evolutionary history is typically found in early-
type (dwarf) galaxies (e.g., Layden & Sarajedini 2000;
Bonifacio et al. 2004; Hendricks et al. 2014), allowing
one to place further constraints on the progenitor’s iden¬
tikit. This emphasises the importance to increase the
density of observational sampling of absorbers within
the inner 30 kpc of the M 31 halo. This will be eased in
the future by the increasing number of known AGN at
suitable redshifts and magnitudes behind the M31-M33
system, as listed, e.g., by Huo et al. (2013).

AK gratefully acknowledges the Deutsche Forschungs-
gemeinschaft for funding from Emmy-Noether grant Ko
4161/1. The authors would like to thank K.R. Sembach
for designing the observing program and D.B. Reitzel for
helpful discussions. We also thank the anonymous referee
for a swift and constructive report.

REFERENCES

Asplund, M., Grevesse, N., & Sauval, A. J. 2005, Cosmic
Abundances as Records of Stellar Evolution and
Nucleosynthesis, 336, 25

Belokurov, V., Evans, N.W., Irwin, M.J., et al. 2007, ApJ, 658,
337

Bernard, E. J., Ferguson, A. M. N., Richardson, J. C., et al. 2015,
MNRAS, 446, 2789

Bonifacio, P., Sbordone, L., Marconi, G., Pasquini, L., & Hill, V.
2004, A&A, 414, 503

Boylan-Kolchin, M., Springel, V., White, S. D. M. & Jenkins, A.
2010, MNRAS, 406, 896

Braun, R., Thilker, D., & Walterbos, R. A. M. 2003, A&A, 406,
829

Braun, R., & Thilker, D. A. 2004, A&A, 417, 421
Brown, T. M., Smith, E., Ferguson, H. C., et al. 2006, ApJ, 652,
323

Danforth, C. W., & Shull, J. M. 2008, ApJ, 679, 194
Danforth, C. W., Keeney, B. A., Stocke, J. T., Shull, J. M., &
Yao, Y. 2010, ApJ, 720, 976

10

A. Koch et al.

Danforth, C.W., Tilton, E.M., Shull, J.M., et al. 2014, ApJ,
submitted (arXiv:1402.2655)

Fardal, M. A., Babul, A., Guhathakurta, P., Gilbert, K. M., &
Dodge, C. 2008, ApJ, 682, L33

Fardal, M. A., Weinberg, M. D., Babul, A., et al. 2013, MNRAS,
434, 2779

Ferguson, A. N. M., Irwin, J. J., Ibata, R. A., Lewis, G. F.,
Tanvir, N. R. 2002, AJ, 124,1452
Gilbert, K. M., Guhathakurta, P., Kollipara, P., et al. 2009, ApJ,
705, 1275

Giroux, M. L. &; Shull, J. M. 1997, AJ, 113, 1505
Grillmair, C. J. 2009, ApJ, 693, 1118

Green, J. C., Froning, C. S., Osterman, S., et al. 2012, ApJ, 744,
60

Haardt, F. & Madau, P. 2012, ApJ, 746, 125
Helmi, A., & White, S. D. M. 1999, MNRAS, 307, 495
Hendricks, B., Koch, A., Walker, M., et al. 2014, A&A, 572, AA82
Huo, Z.-Y., Liu, X.-W., Xiang, M.-S., et al. 2013, AJ, 145, 159
Ibata, R. A., Gilmore, G., Irwin, M. J. 1994, Nature, 370, 194
Ibata, R., Irwin, M., Lewis, G., Ferguson, A. M. N., Tanvir, N.
2001, Nature, 412, 49

Ibata, R., Chapman, S., Ferguson, A. M. N., et al. 2004, MNRAS,
351, 117

Ibata, R., Martin, N. F., Irwin, M., et al. 2007, ApJ, 671, 1591
Ibata, R. A., Lewis, G. F., McConnachie, A. W., et al. 2014, ApJ,
780, 128

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005,
A&cA, 440, 775

Keeney, B. A., Danforth, C. W., Stocke, J. T., Prance, K., &
Green, J. C. 2012, PASP, 124, 830
Keeney, B. A., Stocke, J. T., Rosenberg, J. L., et al. 2013, ApJ,
765, 27

Koch, A., Rich, R. M., Reitzel, D. B., et al. 2008, ApJ, 689, 958
Layden, A. C., & Sarajedini, A. 2000, AJ, 119, 1760
Lehner, N., Howk, J. C., Thom, C., et al. 2012, MNRAS, 424,
2896

Lehner, N., Howk, J. C., & Wakker, B. P. 2015, ApJ, 804, 79
Lewis, G. F., Braun, R., McConnachie, A. W., et al. 2013, ApJ,
763, 4

Mackey, A. D., Huxor, A. P., Ferguson, A. M. N., et al. 2010,

ApJ, 717, Lll

Martin, C., Carlin, J. L., Newberg, H. J. & Grillmair, C. 2013,
ApJ, 765, L39

Mathewson, D. S., Cleary, M. N., Sc Murray, J. D. 1974, ApJ,

190, 291

Mayer, L., Mastropietro, C., Wadsley, J., Stadel, J., Sz Moore, B.
2006, MNRAS, 369, 1021

McConnachie, A. W., Irwin, M. J., Ibata, R. A., et al. 2003,
MNRAS, 343,1335

McConnachie, A. W., Irwin, M.J., Ibata, R.A., et al. 2009,
Nature, 461, 66

Mori, M., Sz Rich, R. M. 2008, ApJ, 674, L77
Mould, J., &; Kristian, J. 1986, ApJ, 305, 591
Newberg, H. J., Yanny, B., &; Willett, B. A., 2009, ApJ, 700, L61
Osterman, S., Green, J., Froning, C., et al. 2011, Ap&SS, 335, 257
Putman, M. E., de Heij, V., Staveley-Smith, L., et al. 2002, AJ,
123, 873

Rao, S. M., Sardane, G., Turnshek, D. A., et al. 2013, MNRAS,
432, 866

Reitzel, D. B., & Guhathakurta, P. 2002, AJ, 124, 234
Rich, R. M., Mighell, K. J., Freedman, W. L., Sz Neill, J. D. 1996,
AJ, 111, 768

Richardson, J. C., Ferguson, A. M. N., Johnson, R. A., et al.

2008, AJ, 135, 1998

Richardson, J. C., Irwin, M.J., McConnachie, A.W., et al. 2011,
ApJ, 732, 76

Richter, P., Westmeier, T., Briins, C. 2005, A&A, 442, L49
Richter, P., Charlton, J. C., Fangano, A. P. M., Bekhti, N. B., &
Masiero, J. R. 2009, ApJ, 695, 1631
Savage, B. D., &; Sembach, K. R. 1991, ApJ, 379, 245
Searle, L. & Zinn, R. 1978, ApJ, 225, 357
Sembach, K. R., &; Savage, B. D. 1992, ApJS, 83, 147
Sembach, K. R., Wakker, B. P., Savage, B. D., et al. 2003, ApJS,
146, 165

Shull, J. M., Stevans, M., Danforth, C., et al. 2011, ApJ, 739, 105
Stocke, J. T., Danforth, C. W., Shull, J. M., Penton, S. V., &
Giroux, M. L. 2007, ApJ, 671, 146
Thilker, D. A., Braun, R., Walterbos, R. A. M., et al. 2004, ApJ,
601, L39

Wakker, B. P., Sz van Woerden, H. 1991, A&A, 250, 509
Wakker, B. P., Sz van Woerden, H. 1997, ARA&A, 35, 217
Yanny, B., et al. 2003, ApJ, 588, 824

York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ,
120, 1579